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Abstract. We derive rigorously explicit formulas of the Casimir free energy 
at finite temperature for massless scalar field and electromagnetic field con- 
fined in a closed rectangular cavity with different boundary conditions by zeta 
regularization method. We study both the low and high temperature expan- 
sions of the free energy. In each case, we write the free energy as a sum of 
a polynomial in temperature plus exponentially decay terms. We show that 
the free energy is always a decreasing function of temperature. In the cases 
of massless scalar field with Dirichlet boundary condition and electromagnetic 
field, the zero temperature Casimir free energy might be positive. In each 
of these cases, there is a unique transition temperature (as a function of the 
side lengths of the cavity) where the Casimir energy change from positive to 
negative. When the space dimension is equal to two and three, we show graph- 
ically the dependence of this transition temperature on the side lengths of the 
cavity. Finally we also show that we can obtain the results for a non-closed 
rectangular cavity by letting the size of some directions of a closed cavity go- 
ing to infinity, and we find that these results agree with the usual integration 
prescription adopted by other authors. 



PACS numbers: ll.lO.Wx 



1. Introduction 



The Casimir effect was predicted in 1948 [T] as an effect due to vacuum fluc- 
tuation of quantum fields. When attempting to calculate the Casimir energy, one 
inevitably faces the problem of summing a divergent series. There have been a 
number of different regularization methods proposed and used to regularize the 
infinite sum to extract a physical finite quantity. Among these methods, zeta reg- 
ularization techniques have been widely used recently. One can see for example 
the articles [3 ll H H Hj and the books by Elizalde et al [i [10] and Kirsten 
[llj . This method has been extended to calculate the Casimir energy at finite 
temperature [12l [IHl [HI [TH Historically, Casimir effect was calculated for 

electromagnetic field confined between two infinitely conducting parallel plates in 
four dimensional space-time. Later on, Casimir energy has been calculated for 
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scalar field, spin 1/2 field and electromagnetic field in more general space-time. 
Among the different geometries of space that have been under consideration, rect- 
angular cavities of different dimensions are among the most extensively studied 
[l3[T2l[li[li[20l[2ll[23[23[2l[25l|2i^ partly due to the simple 

geometry and also the well-developed mathematical tools. Various aspects of the 
effect, such as the low and high temperature expansions of the Casimir energy or 
force [121 nil [m [ITl [SU [32] , the attractive or repulsive nature of the Casimir force 
[m [13 UHl [21 [Ml , the effect of extra dimension [H [TH [T71 [21] , etc, have been 
discussed. 

A p-dimensional rectangular cavity inside a d-dimensional space is a space of 
the form Qp^d = [0, Li] x . . . x [0, Lp] x R^^p. When d = p, we say that the 
cavity is closed, and when p < d, the cavity is non-closed. The paper by Ambj0rn 
and Wolfram [17] can be considered as the pioneer work in the calculation and 
discussion of Casimir effects at finite temperature for massless scalar field and 
electromagnetic field confined within a rectangular cavity. By using dimensional 
regularization technique, they found that the Casimir energy can be expressed 
using Epstein zeta function whose analytic continuation is well-known. Ambj0rn 
and Wolfram were also able to obtain the high and low temperature expansions of 
the free energy by using the Chowla-Selberg formula _34J for Epstein zeta function. 
Their formulas work for p < d, whereas for the case of closed cavity (i.e. the p = d 
case), they modified the p < d formulas to remove divergencies based on physical 
arguments. However, the divergencies for the high and low temperature expansions 
were removed separately and they did not justify that the two results coincide at 
any temperature. 

Special cases of the results of Ambj0rn and Wolfram have been reproduced and 
extended by several authors using zeta regularization or other methods, see for 
examples [l4l[15l[T8l[19l[22l[24l[27l[28l[33]. In particular, there have been an 
extensive study of the Chowla-Selberg formula for general Epstein zeta function [H 
[6l[l2l[l3l[3Tl[3ll[36l[37l[38l[39] with the aim to obtain the low and high temperature 
expansions of the Casimir energy. However, to the best of our knowledge, no one 
has derived the Casimir energy for fields confined in closed rectangular cavities 
correctly (without divergent terms) purely by zeta regularization techniques. One 
can read for example the third paragraph in the introduction of [29], where they 
pointed out this divergency problem in some of the literatures (e.g. [40]). In [29] . 
the authors also mentioned that it is desirable to obtain a closed formula for the 
free energy of the electromagnetic field confined in a three dimensional rectangular 
cavity that is valid for all temperature. 

In this paper, we solve a more general problem. We derive the Casimir free energy 
at finite temperature for massless scalar fields and electromagnetic fields confined in 
a closed rectangular cavity with different boundary conditions, by employing zeta 
regularization techniques. We derive explicit formulas for the free energy, in the low 
and high temperature regions respectively. However, we want to emphasize that 
both the low and high temperature formulas are valid at all temperature. Their 
difference lies in the manifestation of the leading behavior of the free energy at low 
and high temperature respectively. The advantage of using the zeta regularization 
approach is that we can derive formulas that work for any dimension d > 2 at 
one shot. With the further help rendered by the Chowla-Selberg formula, we can 
compute the free energy effectively. We show some results graphically when d — 2 
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and d = 3. On the other hand, we also study some behavior of the free energy 
using the formulas we derive. In particular, we find that the free energy is always 
a decreasing function of temperature. In the cases of massless scalar field with 
periodic and Neumann boundary conditions, the zero temperature free energy is 
always negative. Therefore, the free energy is negative at all temperature. In the 
cases of massless scalar field with Dirichlet boundary condition and electromagnetic 
fields, the zero temperature free energy can be positive. We study the cases when 
d ^ 2 and d = 3, and we leave a more detail study of the general cases to another 
paper. When the zero temperature free energy is positive, we can conclude from the 
decreasing behavior of the free energy that there is a unique transition temperature 
(depending on the side lengths of the cavity) where the sign of the free energy 
changes from positive to negative. We show graphically the dependence of this 
transition temperature on the side lengths when d — 2 and d — 3. In the last section, 
we show how to obtain the corresponding results for a non-closed rectangular cavity 
^p,d by letting the size oi d — p directions of a closed cavity going to infinity. We 
find that our results are in agreement with those based on the method of changing 
the summation in d — p directions to integration, which is commonly adopted by 
other authors. 



2. Casimir Energy at finite temperature 

For a massless scalar field cj) in d-dimcnsional space fl maintained in thermal 
equilibrium at temperature T, the Helmholtz free energy is conventionally defined 
as 

where (3 — l/T and Z is the partition function given by 



(2.1) ^ = nr 



k 

Here Wk is the frequency associated with the eigenmode 0k of the field, and the 
symbol ' in the product means that the term Wk = is to be omitted. More 
precisely, the free energy F is equal to 



(2.2) F = -i log Z = i ^'ck + ^ E' l°g (1 - ^ 



^ 2^ /5 k 



The first term 



2^ 

k 



is the zero temperature contribution to the free energy, also known as Casimir free 
energy. The summation is divergent and regularization is needed to obtain a finite 
value. There are various regularization techniques that have been employed. One 
of the conventional methods is to introduce the zeta function Cn{s) (see e.g. [9|lll)): 

k 

It is well known that Cn{s) can be analytically continued to the complex plane with 
possible simple poles at s = ^5^, I — 0,1,2.... In the case Cn(s) is regular at 
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s = —1/2, we can define 

k ^ ^ 

In general, as was proposed by Blau and Visser [2] , one should introduce a constant 
A with dimension (length)"^ and define 

=i (limi;^ (-l+e) + En { - e 



4 \e^0 

where P.P. means principal part and Eq{s) is the normalized zeta function 



2 N 



X 

k ^ 

Since Cn(s) may have a simple pole at s = —1/2, we can write 
(2.4) Co(.s) = ^+ro + ofs+i 

A straightforward computation gives 

i^o=i(ro + ri logA^). 

If Cn(s) is regular at s = —1/2, ri — ro — (^o(— 1/2) and we get back the 
definition 

The second term in 



k 

is known as the thermal correction to the free energy. Due to the exponential term, 
it is a finite sum. Hence if we are interested in the low temperature behavior of the 
free energy, we can use the expression 

(2.5) = -i logZ = i (ro + ri log A^) + 1 ^'log (l - e-'"'^-) . 

However, this expression is not convenient for studying the high temperature be- 
havior of the free energy. 



> 0. 



Remark 2.1. Differentiate (j2.5p with respect to (3, we find that 



V log (l - e-'3'^M + - V- 



Therefore the free energy is always an increasing function of /3, and thus a decreasing 
function of the temperature T. Hence, if the zero temperature free energy is 
negative, then the free energy F will be negative for all temperature. 



CASIMIR ENERGY IN CLOSED RECTANGULAR CAVITIES 



5 



It has been taken for granted (or taken as definition) that the partition function 
Z can be calculated using the path integral 



(2.6) 



-1/2 



boundary 
conditions 



where 



dt"^ ^ dxf 

is the {d + I)-dimensional Euclidean d-Alembertian operator and /i is a normal- 
ization constant with the dimension of mass. In the imaginary time formalism 
(or Matsubara formalism) of finite temperature field theory, one imposes periodic 
boundary condition with period /? in time direction. In the spatial direction, ip is 
assumed to have the same boundary condition as (f>. The eigenvalues oi -De are 
then given by 

o \ 2 



Using the zeta regularization method, one defines 

(2.7) c(^)= E E' 



O \ 2 

Zirn 



n— — C30 k 



f3 



which is an analytic function of s when Re s > (d -f l)/2. Here the symbol ' in 
the double summation means that a term where (n, uJk) = (0, 0) should be omitted. 
One then analytically continue ({s) to the complex plane and the logarithm of (|2.6p 
is then equal to 

(2.8) logZ = ic'(0) + i(logAi2)C(0). 

Most people set /i^ = 1 and claim that logZ — logZ. However, we are going to 
show that this is not true when there are some modes (/)k with cjk = 0. Since in the 
definition of the partition function (j2.1D . we omit the terms where o^k — 0, therefore 
it is natural to single out the contribution from o^k = terms and write (|2.7p as 

2^\-^^ 



C(s) - 2Af j Cr{2s) + C(s), 

where C-r(s) = S^^^i ^^e Riemann zeta function, Af is the number of modes 

(^k with tJk = and 

C(s) = 




It is well known that C-r(s) has analytic continuation to the whole complex plane 
with single pole at s = I. On the other hand, using standard techniques, for 
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Re s > (d + l)/2, (■(s) is analytic and is given explicitly by 

Here -ftr^(z) is the modified Bessel function of second kind (see e.g., 3.471 in [41j). 
From this, we find that (|2.8p is given by 

(2.9) 



logZ = --log(/3/z)'-/3 ylogM2 + (l-log2)ri + ^ + l ^ log (l 
whereas 

(2.10) logi:=ic'(0) + i(logAi2)C(0) 

= -/3(^logM' + (l-log2)ri + ^ + i ^ log(l 

Compare these expressions with p.Sp . we note that when Af (i.e. in the 
presence of Wk = modes), logZ 7^ logZ, but logZ = logZ if we identify A with 
efi/2. 

It has been noticed by several authors (see e.g. [2l[5l[T4]) that the Casimir energy 
at zero temperature can be defined by 

lim (-ilogZ 

In view of what we have obtained above, due care has to be taken in the presence 
of = modes. In this case, we should replace logZ by logZ. From (|2.5p 
and (|2.10p . we can write 

F = -ilogZ = -i (^logZ + ylog(/3Ai)' 

and therefore 



hm (-- 

P^oc \ (3 



with the identification A = e/i/2. 

The constants /i or A contribute ambiguities to the Casimir free energy. However, 
in most of the cases of interest, the function (n{s) is regular at s = —1/2. This 
is equivalent to ri — 0. Using the zeta function C('S)i we can characterize such 
cases by C(O) = —Af. Hence if C(O) = —A/", the Casimir energy turns out to be 
independent of or A and can be calculated by using 

(2.11) i;^ = _ 1 (^'(o) + 2AAlog/3), 

2/3 

in contrast to the usual prescription F = — ^^^'(0). 

The expression for logZ (|2.9p . with the presence of Wk = terms has been ob- 
tained in [Mj. However, in [M], the discrepancy between Z and the thermodynamic 
partition function Z was not emphasized. On the other hand, a computation sim- 
ilar to what we perform above was done in [3], without taking into consideration 
the Wk = terms. 
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In some of the studies (e.g. ^15,), the (internal) energy E of the system was 
calculated instead of the free energy F. They are related by 

—IP- 

Another important thermodynamic quantity-the entropy S, can be calculated from 
the free energy by the formula 

In view of Remark 12. li it is always non-negative. In the following, we will only 
compute the free energy explicitly. We leave the readers to work out the energy 
and entropy themselves by using these two formulas. 

Remark 2.2. For the sake of convenience of presentation, in this section, we have 
assumed that is a massless scalar field. However, the same reasoning works for 
other quantum fields. 



3. Homogeneous Epstein Zeta function 

Now we want to compute the derivative at zero of the Epstein zeta function using 
Chowla-Selberg formula. In association with the application of zeta regularization 
method, Chowla-Selberg formula has been extensively used to express Epstein zeta 
function in the form which facilitates the study of the function in certain limits 
[llElinilTSlISlllMllMlEllllHlISS]. However, we are unaware of anything done 
regarding the explicit computation of the derivative at zero of the homogeneous 
Epstein zeta function. 

In this paper, we only consider the homogeneous Epstein zeta function in n 
variables in the following form: 



■^_E,n(s; ai, . . . , a„) = ^ 



1 



([aiA;il2 + . . . + lunkJ'^Y ' 
(fci,...,fc„)eZ"\{o} i ij I < I n m J 

This sum is convergent for s > n/2. Under a scaling i-^ Aa^, we have 

(3.1) ZE,n{s; Aai, . . . , Aa„) = A^^*Z£;,„(s; ai, . . . , a„). 

To find the derivative at s = 0, we first derive the Chowla-Selberg formula for 
Epstein zeta function. For fixed 1 < m < n — 1, we can write 

ZE,n{s; ai, . . . , a„) =ZE,7n{s; ai, . . . , «„) 



E E 



([aifcil2 + . . . + [a„fc„12)«' 
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For the second term, we have 
(3.2) 

E 

keZ"x(Z"-'"\{o}) 
1 



([aifci]2 + . . . + [anknYY 



r(s; JO 



ls-\ ^ g-t([ai/cil-+... + [a„fc„]^)g-t([a„+ife„+i]^ + ... + [a„fe„]'')j^ 

keZ'"x(Z"-"\{o}) 



E 



kGZ'"x(Z"-'"\{0}) 



n^liaj r(s 



.TO A -'^ / 

^ E ^n—rn 1 "TT ; ^m+1 : ■ ■ ■ i ) ^" \" n,m('^; ^li ■ ■ ■ i j i 

2 / i (s) 



where 



^n,m ^1 1 ■ ■ • i ^n) 



27r^ 



Em 
3 = 1 



E 

ke(Z'"\{o})x(Z"-"\{o}) 



/ 

[ 1 



E 



n 2 ' 



E [«'^']' 



Combine together, we have the Chowla-Selberg formula 
(3.3) 

ZE.n{s; ai, . . . , a„) = ZE,mis; ai, . . . , am) 

^nr=i«.- r(.) ■ ^ 2 



r(s) 



The function Tn^m{s] oi, . . . , a„) is an analytic function of s on C. Using the fact 
that the Riemann zeta function Cr{s) is meromorphic on C with a single pole 
at s = 1 and the fact that ZE.i{s\a) — 2a^2*Ci?,(2s), we obtain by recursion a 
meromorphic extension of ZE.n{s] oi, . . . , a„) to C with single pole at s = n/2. On 
the other hand, one can also use the Chowla-Selberg formula (|3.3p to prove the 
reflection formula 



(3.4) 



TT T(s)Z£;^„(s;ai, . . . ,a„) 



nn 
3 = 1 



TT' 2 /n \ (n 1 1 ' 

-^l2"'J^^'"l2^^'V---'^, 



by induction (see e.g. [42]). Putting m = 1 and s = in (|3.3p . using the reflection 
formula (^3), the fact that (l/r(s))|^=o = (s/r(s + l))|^=o = and Ci?(0) = -1/2, 
we find that 



(3.5) 



Z£,„(0; ai, . . . , a„) = Ze,i{^\ a) = 2Ci^(0) = -1. 
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On the other hand, if ai < . . . < a„, by putting m = n — 1 in the Chowla-Selberg 
formula (|3.3p . we obtain by recursion 



n-l 



n-l 



1 



T(s) ^ a^«-J 

' j=i "j+i lli=i 



ai 



ai 



kGZJ\{0}ni=l + 



k 



E 



a; 



2iima 



i+i' 



which express the Epstein Zeta function as a sum of Riemann zeta functions plus a 
remainder which is a multi-dimensional series that converges rapidly. This formula 
can be used to effectively compute the Epstein zeta function to any degree of 
accuracy. 

To compute the derivative ^^(s; ai, . . . , a„) at s = 0, we differentiate the 
Chowla-Selberg formula l|3.3p with respect to s and setting s = 0. This gives 

(3.6) 

■^£;,n(0; fli, ■ ■ • , an) 

^^E.mi^' "^l' ■ ■ ■ ' ^ni) ^ VHZ 1 — ZE,n-m (^~"^' ■ ■ ■ , 0,nj + Rn,m{o.l, ■ ■ ■ , a„) 



^^£:.m(0; ai, . . . , am) + 



1 



where 
(3.7) 



Rn.ra (^1 ; ■ ■ ■ i ) 



Em 
.7 = 1 





2 






11^=1 «J 























E 



Afo}) 




Using p.ip and p.Sp . we find that under the scaling i— s- Aa^, we have 
(3.8) Z^„(0;Aai,...,Aa„) = 2 log A + Z^_„(0; ai, . . . , a„). 

4. MASSLESS scalar field inside closed RECTANGULAR CAVITY 

In this section, Casimir energy at finite temperature for massless scalar field 
confined within a closed rectangular cavity of dimension d > 2 will be derived. 
Using the notations in Section [21 the d-dimensional space J7 is the rectangular 
box [0, Li] X ... [0, Ld\ with volume V = Li . . . Ld- Without loss of generality, we 
assume that < Li < . . . < Ld. We are going to consider the following boundary 
conditions for the field 0: A) Periodic boundary condition, B) Dirichlet boundary 
condition, C) Neumann boundary condition. 



A) Periodic Boundary Condition. Consider the periodic boundary condition 



with 4>{xi, . 



,Xd) 



,Xj,. 



,Xd) for all 1 < j < d. In this 
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case, the eigenmodes of (/) are 



0k(x) = e v^^^-^^^;, keZ^ 

The correspondmg zeta function C(s) is 

(m,k)eZ'i+i\{Q} \ \^^ \ W \a/ 

=ZE,d+i s; — , — , ■ ■■ 

\ [j Li Ld 

and there is A/" = 1 zero modes of (f) corresponding to k = 0. By p.Sp . Cp.d(0; Li, . . . ,Ld) 
— 1 = —A/". Therefore by (|2.1ip . the Casimir free energy is given by 

^Hi..-.i.) = 4^««(o;|.^ 

Using (|3.8p . we find that under the simultaneous scahng /3 i-^ A/3, i-^ AL^, the 
free energy Fp{Li, . . . , Ld) transform as 

(4.2) Fp{Li, . . . , Ld) ^ Fp{XLu ■ ■ ■ , XLd) = jFpiLi, . . . , Ld). 

Therefore, when studying the free energy, we can define the scaled variables 

called the scaled temperature and the scaled side lengths of the cavity respectively. 
The function V^^'^Fp{Li, . . . , Ld) is then a function of these scaled variables: 



V'/'^FpiL,, ...,Ld)^ ^i-A^+i ( 0; ^, ^, . . . , ^ ) - i loge 



2tt 2tt 2tt\ 1 

with li . . .Id — ^■ 

The Casimir force on the walls Xj = and Xj ~ Lj is given by 

dF 

(4.3) = -l- 



dL, 



and the corresponding pressure is 
(4.4) P, = 



'~ V ■ 



A.l. Low temperature expansion. By putting to = 1, ai = 27r//3, aj — 
2'K/Lj-i when 2 < j < d+1 in (|3.6p . we obtain the low temperature (T = 1//3 <C 1) 
expansion 

. X . ^ Li...Ld fd+l\ fd+l 

(4.5) Fp(Li, . . . , Lrf) - - ^^iir^r j Z^.d [^-^-.Lu 



E log l-e- 
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which have the form of (|2.5p . We find directly that the zero temperature Casimir 
energy is 

(4.6) FULu . . . , l.) = - ^^^-^r z,., [^■,L„ ...,1 

which agrees with (3.4) in [17]. A similar result was obtained by Edery [27] using 
multidimensional cut-off technique. By the definition of the Epstein zeta function, 
the term ()4.6|) is strictly negative. Remark |2. II then implies the Casimir free energy 
is then always negative for all temperature. On the other hand, we can compute 
an explicit upper bound for the thermal correction term: 



|AFp(Li,...,Lrf)| = 



kez<i\{o} 



log 1 - e 



+...+ 



1 

< - 



1 — e ^-^d 



d+i ' 



which is an exponentially decay term as /3 — > cxd. 

From (|3.ip and ()4.6|) (or by (14. 2p ). we see that under the space scaling Li 
XLi, 1 < X < d, the zero temperature free energy transforms as 



(4.7) 



Namely, the zero temperature free energy is inversely proportional to the dimension 
of space. This scaling property breaks down at positive temperature. However, (|4.2p 
shows that this scaling behavior will hold if the temperature is also scaled inversely. 
On the other hand, differentiating the equation on the right hand side of (|4.7p with 
respect to A and setting A = 1, we get 



dF" 



= -F" 



From the definition of pressure (j4.4p . we find that at zero temperature, the equation 
of state 



(4.8) 



F^ 



(Pi 



+ Pd)V 



holds. When the cavity is a hypercube (i.e. when Li = . . . = Lj), this implies 
that the zero temperature free energy F^ always has the same sign as the force 
and pressure. At finite temperature, as a correction to (|4.8p . (|4.2p gives us the 
well-known thermodynamic relation 



(4.9) 



F 



, dF dF OF , 

oLi oLd op 



Pd)V - TS, 



where S is the entropy (|2.13p . 

Using Arithmetic-Geometric inequality, we find that when V = Li 



. L(i is fixed. 



Li 



+ {LdkdY > dih ■ 
^V>d(fci.. 

l^d J 



..kd)'iV'i, 



and equalities hold if and only if Li = L2 = ■ ■ ■ — Ld- Therefore, we conclude 
from (|4.5p that at fixed volume, the Casimir energy achieved its maximum when 
Li = L2 — ■ ■ ■ — Ld- 
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A. 2. High temperature expansion. By putting m = d, aj = 27r/ij, I < j < d, 
ad+i = 27r//3 in p.6p . we obtain the high temperature (T = 3> 1) expansion of 
the free enrgy 



(4.10) 



d + l / , , \ 
TT — ^ .„/^d+l\ 1 



Fp(ii, . . . , Ld) = -^^ii • • • j Cfl(d + 1) - ^ log(2^/3) 




27rTO 



2 

J J 



The leading term 

is the usual Stefan-Boltzmann term. In some of the existing literature (e.g. [2]), 
the second leading term i log(27r/3) was overlooked. However, since this term does 
not depend on the dimension of the space Li, . . . , Ld, it does not contribute to the 
Casimir force. Nevertheless, this term is essential for the validity of the thermody- 
namic relation (|4.9p . The last term in (|4.10p is an exponentially decay term. More 
precisely, it is bounded above by 



2-kLi 
1 -|- e~ fiVd 



Ambj0rn and Wolfram obtained a similar high temperature expansion in |17| (see 
(7.10)). They considered the non-closed cavity case and let p = d in the formula 
valid for p < d, and then removed the divergent term by subtracting the free bose 
gas result. They did not justify their result mathematically. Here we have proved 
this formula rigorously. 

We would also like to mention that the general structure of the high temperature 
expansion of free energy of gases inside cavities in curved space-time has been 
calculated (see e.g. [131 1331 US]). Our result here can be considered as special case 
of their result. 

B) Dirichlet and Neumann Boundary Condition 

B.l. Dirichlet Boundary Condition. The eigenmodes of satisfying the Dirich- 
let boundary condition 0(x)|gj2 are 

0k(x) = nsin('^x,V keN"^. 
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The corresponding zeta function C(s) is 



CD,d{s; Li, . . . , Ld) — 

(m,k)eZxN'' 



P J \LiJ \L, 
There is no zero mode of 6 in this case. 



B.2. Neumann Boundary Condition. For the Neumann boundary condition 
(9n</'(x)|gQ — 0; where n denotes the unit vector normal to the surface dQ, the 
eigenmodes of (j) are 

M^)^llcos(^xX ke(NU{0})^ 
The corresponding zeta function C(s) is 

CAf,d(s;ilj • ■ • j-^d) = ^ ' 

(m,k)eZx(NU{0})'' 

There is A/" = 1 zero mode of (f> in this case corresponding to k = 0. 
Since 

9{ki, ■■■,kd) ^2-" 5] (1 - <5fei,o) ... (1 - -Jfe^o) 5(^1, . . . , fed), 

^ 5(A;i,...,fcrf)=2-'' ^ (l + 4i.o)...(l + 4,.o)5(fci,...,M 

ke(Nu{o})'' keZ'i 

for any function g satisfying g(fci, . . . , — fc^, . . . , /c^) = 5(^1, ■ ■ ■ , ki, . . . , kd), 1 < i < 
d, we have 

Cd/nA-'^-^ Li,---, Ld) - 2-" (^{Tlf ) Cfl(2s) 



^ — 1 l<mi < . . .<mj <(i 



From this, it is easy to check that CD.d(0; Li, . . . , Ld) — and CAf,d(0; Li, . . . , Ld) 
— 1. Therefore, by (|2.1ip the free energy is given by 



(4.12) 



27r TT 



J = l l<mi<...<mj<d ^ ^ ^ 

+ ^ log/3 -^^/A'^log/3> 

where = and Qjq — \. Compare to the free energy of the periodic case (|4.ip . 
we have 

d 

(4.13) Fc/jv(ii,...,^d)=2-'^E(Tl)'"-'' E Fp(2L™,,...,2L™J. 

J — 1 1 < ?7i 1 < . . . < m J < d 
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Figure 1: The graph on the left shows the free energy Fd{Li, L2) as a function of Li 
when V = L1L2 = 1, at T = 0, 0.3, 0.6, 0.9, 1.2, 1.5. The graph on the right shows the free 
energy Fd{Li,L2) as a function of T when Li = 0.4, 0.55, 0.7, 0.85, 1.0 and V = L1L2 = 1. 




Using this formula and (|4.2p . we find that under the simultaneous space-time scaling 
/3 ^—^ A/3, Li I— > XLi ,1 < i < d, the free energy for the Dirichlet and Neumann 
conditions Fiy/j^{Li, . . . ,Ld) behave in the same way as the free energy for the 
periodic condition Fp{Li, . . . , Ld) (14. 2p . and thus the thermodynamic relation (|4.9|) 
also holds in these cases. 

The low and high temperature expansions of the Casimir free energy Fjj/jy can 
be obtained directly from (|4.13p and the corresponding expansion for Fp. As in 
the periodic case, the zero temperature free energy for the Neumann case is always 
negative. However, the sign of the zero temperature free energy of the Dirichlet 
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Figure 3: The free energy Fd{Li, L2, L3) as a function of T when V — L1L2L3 — 1 and 
k2 — L2/L1 = 1, 5 at fca = I/3/I/1 = 1, 5, 10, 20 respectively. 




Figure 4: The free energy Fd{Li, L2, L3) as a function of fca — L3/L1 when T — 0,2,3 
and V = L1L2L3 = 1 at = L2/L1 = 1,2,4,8. 



case depends on the parameters Li, . . . ,Ld. There have been a lots of discussions 
about this in the literature, see e.g. [TSl [T71 [TFl ■ By Remark 1 2. 11 we know that 
fixing Li, . . . , Ld, if is positive, then there exists a unique T — T{Li, . . . , Ld) 
such that Fd(Li, . . . , Ld) change from negative to positive. The scaling property 
of free energy (g^ shows that T{XLi, . . . , XLd) = \^^T{Li, ...,Ld). We study 
this transition temperature graphically for d — 2 and d = 3 (see Figure 5, 6). 
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Figure 5: The transition temperature T{L\,L2) for Fd{L\, L2) as a function of Li when 
V = L1L2 = 1. F^{Li,L2) is positive when 0.6045^ < L1/I/2 < 1/0.6045^. 




0.6 0.8 1 1.2 1.4 1.6 1.8 4 6 8 10 12 14 16 18 20 

k2 k3 



Figure 6: Left: When 0.5733 < k2 = L2/L1 < 1.7444, there is a unique fca = ^3(^2) 
such that Fd{Li, L2, Ls) > for all L3/L1 > fca. The graph shows logfcs as a function of 
k2. Right: The transition temperature T{Li, L2, L3) for Fd{Li, L2, L3) as a function of 
ks = L3/L1 when V = 1 and fo = L2/L1 = 0.75, 1, 1.25, 1.5. 



Table 1 



^2 — L2/L1 


The range of ^3 — L3/L1 where F'^{Li, L2, L3) > 


0.75 


ks > 4.4972 


1.0 


ks > 4.2471 


1.25 


ks > 5.2999 


1.5 


k3 > 8.8571 



In the high temperature regime, the leading term is 



-^^d+i 



It comes from the j = d term in (|4.13p and it is again the Stefan-Boltzman term as 
in the periodic case (14. lip . The term proportional to ^ log (3 is also present but in 
the Dirichlet case, its sign depends on d. Unlike the periodic boundary case (|4.10p . 
now we have terms proportional to 1//3-' for all 1 < j < d + 1. 
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5. Massless vector field (electromagnetic field) 

As discussed in [17 , for massless vector (spin 1) field (or electromagnetic field) 
inside a d-dimensional space $7, the field strength is represented by a totally anti- 
symmetric rank-2 tensor F**" satisfying the equations 



=0, 



where FMi - Md-i = ^tii,...,tJ.d-i,i^,>^p^^ jg ^jjg ^^^^i tensor of i^'^" and j'^ is the current. 
In the vacuum state, = 0. 

C.l. Perfectly Conducting Walls. In the case that — [0, Li] x . . . x [0, Ld] is a 

rectangular cavity with walls of infinite conductivity, the field satisfies the boundary 
condition 



= 0, 



where is the unit vector normal to the walls dfl and hq = 0. Introducing the 
potentials so that 

F^"" = d'^A" - a" = do,d' = -9,, l<i<d 

and working in the radiation gauge with gauge condition 

(5.1) A° =0, 

^^A' =0, 

we find that the modes of the potentials are given by 

d 



Al 



-ai cos 



1 < i < d. 



where k e (N U {0}) 

The gauge condition (|5.ip implies 

(5-2) E'' 

1=1 



d 



It is easy to see that if two of the /c^'s are zero, Ak = (A^, . . . , A^) is identically 
0. On the other hand, if only a single hi is zero, then for j ^ i, AL = and 



is trivially satisfied. When all fc^'s are nonzero, (|5.2p implies that there is 
a (d — 1) degree of freedom for the vector a — {ai, . . . , ad) for any fixed k S N"*. 
Therefore the zeta function for electromagnetic field confined in rectangular cavities 
with perfectly conducting walls is related to the zeta function for massless scalar 
field under Dirichlct boundary condition by 



CAc,<i(s;ii, • ■ ■ ,^d) =(rf- l)CD,d(s;Li, . . . ,_Lrf) 

, CD,d-i{s; Li, . . . , Lj_i, Lj+i, . 



d 
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There is no ujk = mode and CAc,d(0;ii, . . . ,Ld) — 0. The corresponding free 
energy is 

(5.3) 

d 

Fac (^1; ■ ■ ■ I Fd) — {d~ 1)Fd(Li, . . . , Ld) + ^ Fd{Li, . . . , Lj_i, L^+i, . . . , Ld) 

i=i 

d 

j—l l<mi<...<rnj<d 

C.2. Infinitely Permeable Walls. In the case that fl — [0, Li] x . . . x [0, Ld] is 

a rectangular cavity where the walls are infinitely permeable, the field satisfies the 
boundary condition 

Working in the radiation gauge (|5.ip , the modes of the potentials are given by 
Al =7. sin (^^x.) n cos (^^x,^ e"- 1 < ^ < d, 



d / , X 2 



where k e (N U {0})'', 

The gauge condition (|5.f |) implies 

(5.4) E ^ = 0. 

■i— 1 * 

If all the fci's are zero, Ak = (A^, . . . , yl^) is identically 0. On the other hand, for 
1 ^ j < fixing 1 < ri < . . . < rd-j < rf, let 1 < toi < . . . < < d be such 
that {mi, . . . , mj,ri, . . . , rrf_j} = {f , 2, . . . , d}. If fc^ = • ■ . kr^_- = and km^ ^ 

0,...,k,n, ^0, then = for 1 < Z < d- j and dE!]) reduces to XiLi ^I'^"' = 0. 
This implies that there is a {j — 1) degrees of freedom for the vector {jmi , ■ • • , 7mj ) 
for any fixed {kmi , ■ ■ ■ , km^ ) €W ■ Therefore the zeta function for electromagnetic 
field confined in a closed rectangular cavity with infinitely permeable walls is related 
to the zeta function for massless scalar field under Dirichlet boundary condition by 

d 

CAB,d(5;^l, • • ■ ,^<i) =E(-^' ~ E Cdj (s;imi, . . . ,imj) . 

J— 2 l<mi <. . . <mj <d 

There is no Wk = mode and C,AB,d{Q'T Li, . . . , Ld) — 0. The corresponding free 
energy is (see the detail computation in the appendix): 

d 

(5.5) FAALi,---,Ld}=Y.(j~l) E FD{s;Lrn„...,Lrn^) . 

j—2 l<mi<...<7nj<d 
d 

=2-'*^(2j-d-f) Fp(2L„,,...,2L™J. 

j — l l<7ni <. . . <mj <(i 
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Figure 7: The free energy Faq{Li, L2, L3) as a function of T when V — L1L2L3 = 1 and 
k2 = L2/L1 = 1, 5 at fes = L3/L1 = 1, 5, 10, 20 respectively. 



Notice that when d = 2, Fac{LuL2) = Fn{Li,L2), Fas{Li,L2) = Fd{Li,L2) 
and when d = 3, Fac{Li, L2, L3) = (^i, -^2, -^3)- 

Under the simultaneous space-time scahng (3 1-^ A/3, 1— > XLi, I < i < d, both 
Fac (-^1: • ■ ■ : -^d) and Fab {Li, . . . , Ld) transform as 

FAc^ni^l, ■■■,Ld)>-^ X^^Fac/b{Li, ■ ■ .,Ld), 

and thus the thermodynamic relation (|4.9p holds. 

The low and high temperature expansions of the free energy Fa^/^ can be ob- 
tained from the corresponding expansion of Fp using (|5.3p and (|5.5p . The sign of 
the zero temperature energy F'j^^^^ (ii, . . . , Ld) also depends on the relative size of 
Li, . . . , Ld- 



Table 2 


fc2 = i2/il 


The range of = L3/L1 where L2, L3) > 


0.75 


0.3555 <k3< 2.7033 


1. 


0.4083 < /ca < 3.4298 


1.25 


0.4580 < /c3 < 3.6219 


1.5 


0.5057 < /c3 < 3.4957 



In the high temperature regime, the leading term is 

-(^-i)^r(i±i)c.«i-M). 

which is {d — 1) times the leading term in the scalar field case. This is due to the 
fact that electromagnetic field in d+1 dimensional space-time has d—1 polarization 
states. The term proportional to ^log/3 still present. When d is even, there are 
terms proportional to for all 1 < J < + 1. When d is odd, there are terms 
proportional to 4j for all 1 < j < d + 1 except for j = ^i^. 
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Figure 9: The free energy Fac{Li, L2, L3) as a function of fcs — L3/L1 when T — 0,2,4 
and V = LiL2L:i = 1 at = L2/L1 = 1,2,4,8. 
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Figure 10: Left: When 0.2761 < fc2 = L2/L1 < 3.6195, there exist r?i(fe2) and ??2(ifc2) such 
that Fac{L\, L2, L3) > for all r;i(fc2) < L3/L1 < ri2(k2). The graph shows r]i and 772 as 
functions of ^2. Right: The transition temperature T[L\, L2, L^) for Fac{Li, L2, L3) as 
a function of fcg — L3/L1 when V = 1 and = 1, 1.5, 2. 



When d = 3, we find that the zero point energy i^^^ (Li, L2, -^^3) is given by 



L1L2L3 



167r2 



y ^ 



1 



1 



1 



48 y Li L2 L2 



which is a well known result (see e.g., [23 )■ We show graphically some particular 
values of the transition temperature T{Li, L2, Ld) for Faq {Li, L2, L^) in Figure 10. 
On the other hand, the high temperature expansion of Fac (-^i: ^2, -^3) = Fab {Li, ^2: 

is 



(5.6) 



Fac{Li, L2, L3 
1 



Fab {Li,L2, L3) 



TT^ L1L2L3 TT Ll + L2 + L3, 



213 " ^^^'3 (0; L^\L2\L^') - ^ log(87rLiL2i3) 

LiL^Lz 
2/32 

^1^2-^3 



E 



45 /34 



12 



1 



1 °° 1 
— Y- 



fc=l 



1 



1 



1 



e " 



1 



e -3 



1 



e f> 



1 



Our result gives the correct high temperature limit stipulated by Ambj0rn and 
Wolfram [T7j (equation (7.12)). However, they only obtained the first three terms. 
To the best of our knowledge, we are not aware of any existing study that calculate 
the high temperature limit to the degree of accuracy obtained here. We would 
like to emphasize that the formula (|5.6p is valid for all temperature. In [25], the 
authors calculate this free energy by a different method. They gave the same first 
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two leading terms as above, and no explicit formulas for the remaining terms are 
given. 



There exist many papers on the Casimir energy of massless scalar field or Casimir 
energy of electromagnetic field confined in a p-dimensional rectangular cavity in a 
d-dimensional space [HllIllIIlllIllIiniEniEaEllESlEelE^ p-dimensional 
rectangular cavity in a d-dimensional space is a space of the form flp^^ = [0, Li] x 
. . . X [0, Lp] X E**^^, where < p < d. It can be considered as the limiting case 
of the closed cavity where Li, . . . , Lp <C ip+i = ■■■ = Ld = L or Lj oo for 
p + I < j < d. In the existing literature, when calculating the Casimir free energy, 
usually after setting up the zeta function over in a suitable set, the 

summation over kp^i, . . . , kd is changed to integration. From the mathematical 
point of view, this is not a rigorous treatment since the summation expression 
for the zeta function only works for Re s > |, which does not include the point 
s — 0. To justify this procedure, one actually need to justify that the processes 
of taking analytic continuation and taking limit Lj —^ oo for p + 1 < j < d 
can be interchanged. In this section, we will directly take the limit Lj — > oo for 
p+1 < j < din the expression of Casimir energy for fields inside a closed rectangular 
cavity to obtain the energy of the fields inside a non-closed rectangular cavity. To 
be more precise, the limit when Lj — > oo for p + 1 < j < d of the free energy 
F{Li, . . . , Ld) is always infinite. Therefore we shall consider the free energy density 
/ defined as the limit 



In the following, we assume that < p < d ~ I. By putting m — d — p, aj — 
2TT/Lp+j, 1 < j < d—p, Qd-p+i = 27r//3 and aj — 2'K/Lj-d+p-i,d—p+2 < j < d+1 
in (|3.6p . we find that the free energy Fp{Li, . . . , Ld) (|4.ip is equal to 



where . . . , a„) is defined by p.7p . Now the last term goes to zero as 

Lj~>oofoTp+l<j<d (see appendix). Therefore, 



6. From closed cavity to general case 



fd{Li, . . . , Lp) — 



lim 

Li — >co 
p+l<i<d 





(6.1) 



fp,diLi, . . . , Lp) 




'P+i 



1 



, Ld 



E,d~p 



0; 




27r 27r 



) 
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Next we want to show that the first term in ()6.ip is also zero, i.e. we need to show 
that 



Hm ai . . . a^Z'^ „ (0; ai, . . . , a„) = 0. 

ai—>-0 ' 
l<i<n 



When n = 1, we have 



Z' i(0;a) = 21og 



therefore hnia^o [o^b i(0;a)] — 0. When n > 1, equation (|3.6p gives 



— Zb_„_i (0; ai, . . . , a„_i) + 



nn— 1 
i=l 



■Ci?("-) + Rn,n-i{ai, ■ ■ • ,a„). 



Using a similar argument as before (see appendix), 



lim fli . . . a„„ii?„.„_i(ai, . . . ,a„) = 0. 

l<i<ri-l 



On the other hand, it is obvious that 



hm fli . . . a„ 

l<i<n-l 



nn — 1 
.7 = 1 «J 



Ci?,(n) - 0. 



Therefore, 



lim ai . . . a„Z^ „ (0; ai, . . . , a„) = lim ai . . . a„2'^ „._i (0; oi, 

l<i<n l<i<n 



and we obtain by induction on n that this is zero. Consequently, we find fi'om (|6.ip 
that 



(6.2) /p,d(ii, . . . , ip) = 3Tl~^i ■ ■ ■ LpZE,p+i 

2tt 2 



d+ 1 



; /3, Li, . . . , Lj 



and this agrees with the result in [17j obtained by dimensional regularization 
method. Notice that the right hand side of (|6.2p is not defined when p = d. 
Under the simultaneous space-time scaling (3 A/3, Li \Li ,1 < i < d, the free 
energy density fp,d{,Li, . . . , Lp) transforms as 

(6.3) IpALi. • • • , ip) ^ y-'^fpALi. • ■ • , Lp)- 

Now using the fact that 

Fp {2Lmi , . . . , 2L,nj , 2Lp+i , . . . , 2id) 



lim 

Lj — >co 
p+l<j<d 

^2'^~P lim 

Lj — >oo 
p+l<j<d 



Lp+1 ■ ■ ■ Ld 



Fp ylLyyii J • • • J 2Lm . , 2L„+i , 



,2L, 



(2Lp+i) . . . {2Ld) 



2 ^/p,d (2imi , . . . , 2L„ 



and (|4.13p . (|5.3p . we find that the free energy densities fD,d{Li, . . . , Lp), fN,d{Fi, ■ ■ ■ , Lp), 
fAc,diLi, ■ ■ • , Lp), fAB,d{Li, . . • , Lp) for massless scalar field under Dirichlet and 
Neumann boundary conditions and for electromagnetic field confined in cavity with 
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perfectly conducting walls and with infinitely permeable walls are related to the free 
energy for massless scalar field under periodic condition by: 

(6.4) 

p 

fD/N,d{Ll,---,Lp)—2 ^^(T1)^ ^ ^ fp,j+d-p{'^Ljn^, . . . ,2Lra.) 

j=0 l<rni<...<rnj<p 

/-4c/B.d(-^li ■ ■ ■ ' ^p) 
P 

j—O l<mi<...<mj<p 

The scaling behavior of the free energy density in these cases is the same as the 
periodic case (I6.3p . 

When p = 0, we obtain the vacuum energy of free massless scalar field and 
electromagnetic field in M.'^: 

(6.5) /P,d = fD,d = fN.d = - ^^^\]+^ U{d + 1), 

r i^) 1 

/ac/b,. =-(^-1) l^ ^d+M d+l), 

which are the Stefan-Boltzmann terms. These equations are reasonable since when 
extending to the space R"*, the boundary disappears and the vacuum energy should 
be the same no matter what boundary conditions we start with. 

6.1. Low Temperature Expansion. 

When 1 < p < d — 1, by putting m — p, ai = Li, ^ ^ i ^ p, flp+i = /3 in the 
Chowla-Selberg formula p.3p . we obtain the low temperature (T <C 1) expansion 
of the free energy density (|6.2p : 

(6.6) 



fp,d{Li, . . . , Lp) 



27r— 
2 



(^) 



(jd-p+l 



E 



"1=1 kezp\{0} 



TO 2 \ . , 



27rm/3 



The first term gives the zero temperature energy density and the sum of the last 
two terms is the thermal correction. Note that now the thermal correction contains 
a term proportional to /3 2 . As usual the last term decays exponentially. We 
show in the appendix that the sum of the thermal correction is equal to 



1 



log 1 



dwY . . . dwd-p, 



in agreement with the usual integration prescription to obtain the limit Lj — > 00 
for p + I < j < d. From this formula, we can verify as in the closed cavity case 
that the free energy density is a decreasing function of temperature. On the other 
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hand, ()6.2|) implies that the Casimir free energy is negative at aU temperature for 
all p and d such that < p < d. 

Compare to (j4.5p . we find that we cannot simply set p = d in (|6.6p to obtain 
the free energy in the closed cavity case (I4.5P due to the second term. In fact, by 
using physics argument, Ambj0rn and Wolfram |17| has argued that in order to 
obtain the free energy for closed cavity from this formula, it is necessary to omit 
the second term. 

Using (|6.4p and (|6.6p . one can also obtain the low temperature expansion of 
the free energy densities fo/N.d and fAc/B ^ ior 1 < p < d — 1. We find that in 
the case of scalar field with Dirichlet boundary condition, the thermal correction 
is an exponentially decay term, whereas for the scalar field with Neumann bound- 
ary condition and also for electromagnetic field confined in a cavity with infinitely 
permeable walls, there is an extra term proportional to /3 5 and for the elec- 
tromagnetic field confined in a cavity with perfectly conducting walls, this extra 
term only present when p = I. Just like the periodic case, we can show that the 
thermal corrections are equal to 



(6.7) 



Y: MHc(k) log 1 - e-'^V^?- L-fJ +l-n . . . dw,.,, 



which is in agreement with the usual integration prescription. Here BC = D, iV, Ac, Ab 
and 



Mjv(k) 





fl, 


if k e W, 


-\ 


10, 


otherwise, 


= 1 


Vke (NU{0}f, 






- 1, if fcj 7^ for aU 1 < i < p. 


■1 


I'" 


if fci = for some i, and kj 






otherwise. 



Mab (k) — d — p + j — 1, if exactly j of the ki, . . . ,kp are nonzero. 



From this, we can also conclude that the free energy density is a decreasing function 
of temperature. In the case of scalar field with Neumann condition, we can even 
generalize the conclusion to that the Casimir energy is always negative. However, 
in the case of scalar field with Dirichlet condition and the cases of electromagnetic 
fields, the sign of the Casimir free energy depends on p,d,T and the values of 
Li, . . . , Lp. There have been some discussions on this point in [TSl [ISl [Ml I33j . 

6.2. High Temperature Expansion. 

When 1 < p < d — 1, hy putting m — 1, ai = f3, ai — 2 < i < p + 1, in the 

Chowla-Selberg formula p.3p . wc obtain the high temperature (T ^ 1) expansion 
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of the free energy density ()6.2p : 

(D-«j JpA^i^ ...,Lp)- ■ ■ ■ — — 

r(f) 1 fd 

27r2 p \^ 

2Li . . . d I rr ; 12 \ / 27rm 



P kezp\{0}m=i \i=i / \ \ j=i 



E t^^-^. 



2 



which agrees with the result obtained in [T7]. The leading term is the Stefan- 
Boltzmann term which is equal to the vacuum energy of M.'^ (j6.5p . The second term 
is of order and it is divergent for p ~ d. In [17| . Ambj0rn and Wolfram argued 
that to obtain the p — d case from this formula, one needs to remove the divergence 
by subtracting the free Bose gas result, i.e. replace the second term by 




) ^r^FTT — Ze.i ( ;/3 

TT 



1 



- - ^ {Z'e.p (0; . . . , L-^) - Z'e,,{0, r')) ■ 

Comparing to (I4.10p . we have shown mathematically that this is indeed the case. 
Using (|6.4p and (|6.8p . one can also obtain the high temperature expansion of the 
free energy densities fD/N,d and fAc/B,d for 1 < p < d— 1. We find that the leading 
term for all the cases is equal to the vacuum energy of M.'^ (|6.5p . In the cases of 
Dirichlet and Neumann conditions, there are terms proportional to (3^^ for every 
d — p+l<j<d+la,s well as for j — I. For electromagnetic field, when d is odd 
and p > (d + l)/2, there is no term proportional to /3 ~ . 

7. Conclusion 

We have provided a rigorous derivation of the Casimir free energy at finite tem- 
perature for massless scalar fields and electromagnetic field confined in a closed rect- 
angular cavity with different boundary conditions by zeta regularization method. 
By applying Chowla-Selberg formula, we obtained explicit formulas for the low and 
high temperature expansions of the free energy, which can be written as a sum 
of polynomial order terms in T or plus an exponentially decay term. To the 
best of our knowledge, such explicit formulas for the low and high temperature 
expansions of the free energy of fields confined within closed cavities has not been 
obtained previously. 

We noted that for all the cases considered, the free energy at finite temperature 
F(P; Li, . . . , Ld) transforms as 

F(/3; Li, . . . , Ld) ^ F{\p- Aii, . . . , XL^) - \-^F{(3- ii, . . . , L^), 

under the simultaneous space-time scaling [3 ^ A/3, Li i— > XLi, 1 < i < d. This in 
turn implies the thermodynamic relation 

F={Pi + ...+ Pd)V - TS, 

which has not been observed. 
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On the other hand, we also show that the free energy in aU the cases considered 
is a decreasing function of temperature. For massless scalar field under periodic and 
Neumann boundary conditions, the free energy is negative for all temperature. For 
massless scalar field under Dirichlet boundary condition and for electromagnetic 
fields, the free energy might be positive at zero temperature. When this happens, 
there is a unique transition temperature at which the free energy change from 
positive to negative. This transition temperature is shown graphically for d = 2 
and d — 3. We believe that for massless scalar field under Dirichlet boundary 
condition and for electromagnetic fields, when d > 4, the zero temperature free 
energy will also be positive for (Li, . . . , L^) lying in some domain of R'^. A detail 
study of this is left to another paper. 

In the last section, we show how the free energy for a non-closed rectangular 
cavity can be obtained by letting the size of some directions of a closed cavity going 
to infinity. We prove that the results are in agreement with that based on the 
integration prescription usually adopted by other authors. 

We remark that the discussion given in this paper focused mainly on the low 
and high temperature expansions of the free energy and the properties of the free 
energy. We have not dealt with other thermodynamic quantities such as the force, 
pressure, internal energy and entropy. We hope to consider these quantities in a 
future work. 

Finally, we would like to point out that there exist some controversies regarding 
imposing boundary conditions on a quantum field. Deutsch and Candelas ^46J were 
the first to study the nonintegrable divergences in the renormalized energy density 
near boundaries. This problem has been re-examined by Baacke and Kriisemann 
[¥7] and analyzed in detail recently by Jaffe [IHIIH] and Graham et al |^50, 51, 52 , 5_3]. 
These authors showed that the imposition of boundary conditions on quantum 
fields in Casimir effect calculations leads to non-renormalizable infinities. As a 
result, fixing boundary conditions ab initio invariably results in divergences which 
cannot be removed by renormalization. Basically this problem for electromagnetic 
field with Dirichlet boundary condition can be stated as that no real material is 
perfectly conducting at arbitrary high frequencies. In order to overcome this serious 
problem, Graham and collaborators have developed a new approach which replaces 
the boundary condition by a rcnormalizable coupling between the fluctuating field 
and a non-dynamical background field representing the material. On the other 
hand, there were responses from Milton ,54j. Fulling [Snj and Elizalde [5^ with 
various attempts to resolve this issue. Here we would like to mention the effort 
by Elizalde who has tried to explain the presence of infinities as a result of drastic 
reduction of eigenstates when boundary condition is imposed. He has proposed to 
complement the zeta function method with the Hadamard regularization in order to 
make sense of infinities present in the boundary value problems in Casimir energy 
calculations. However such an approach cannot be taken as a substitute of the more 
physical treatment given in ref . [48l |49l [50l [5TJ [52l [53] . The system considered in 
this paper can be regarded as ideal cases, for which zeta function technique is 
still a useful tool for regularization of vacuum energy density. For a more physical 
treatment, one has no choice but have to take into account of the problem of singular 
behavior near a boundary. 
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Appendix A. 

In this appendix, we gather some mathematical formulas and estimates that we 
need. 

1. We want to prove (|5.3p . By equation (|4.13|) . we find that 

d 

FAb (Li, . . . , Ld) ^'^Cj-d ^ Fp [L,ni, ■ ■ ■ , Lm^) , 

j—1 l<mi<...<mj<d 



where 
Now we compute Cj-d- 

fc=0 ^ ^ 

\fc=o 



k—j 



d-j 
k 



2-' + (j-i)E(-i)' 

fe=0 



d-J 
k 



i-fe 



(d-j) 



id-j) 



k=l 

d-j-1 



fc=0 



d-3 
k 



fc=0 

1 

'--2 



d-J -I 
k 



d-J 



fc=0 



k id- j 
k 



d-]-l 



+ (j-l)(l-2 



=2-\2j-d-l). 
2. We want to show that 



lim 

ai-tO 
l<z<m 



J = l 



Rn,m{0'i, . . . , an) — 0, 



with Rn,m defined by p.7p . Without loss of generality, we assume that ai < . . . < 



a„. Define ai{k) = \/l]JLi ^ , a2(k) = -\/Z]"=m+i K^jP- Then by (P?7)) and 
using 
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where 



we have 



n 



"2j-l" 



ke(z-\{o})x(Z"-'"\{o}) 



Using the inequahty 



n 1 / " \ 



we have 



E 



k6(Z™\{0})x(Z"-™\{0}) 

< E 

kG(Z™\{0})x(Z" — \{0}) 



E 



(^m+l ) • ■■i^n 



<2mam ^ 



1 -|- g a^n-/^ 



a2(k)- 

-\{0} 

m— 1 



2e' 



_ 2x0200 

2 g "^ITl 



E 



a2(k) + 2 e ""•^ 



1 _ e omv^ (fe„.+i,...,fe„)eZ"-'"\{o} 



Prom this, it is easily seen that as — > for 1 < flj < m, Rn,m ~* 0. 
3. We want to show that the integral 



' [ Vlogfl-e-'^^/^MW 



+ |wP 



dwi . . . dwd-p 



is equal to 
(A.2) 

( d-p+l \ 



TT 2 



TO=i kezp\{o} ^ \j=i 
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We split / into two terms Ii and where /i corresponds to k = term and I2 
contains the k e Z*' \ {0} terms. We have 



' ^-p-' y - — dw 



r{d-p) 



V ^ 



Using the formula r(2z) = 22^-i7r-i/2r(z)r (z + i) (8.335 of [H]), we find that h 
is equal to 

p / d-p+i \ 

T^^—pd-p+l 



For I2, set u(k) = yYfj^i [^^] ' "^^ have 



kGZP\{0} 



r (^^) (27r)<'-P/3 ^^^^^j Jo 

Now using the substitution u = Vv^ + w^ and the formula 4 of 3.389 in [41], we 
have 



d-p-i _i 1 d-p+i f d~ p\ 

=2 ^ TT 2 ^^——y 2 r — — {mpv) . 

(mp) 2 \ z / 2 

Combining together we find that I2 is equal to the second term in (|A.2p . thus 
proving our claim. 
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